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In the light of recent interesting experimental work on MgV2 04 we employ the density functional 
theory to investigate the crucial role played by different interaction parameters in deciding its elec- 
tronic and magnetic properties. The strong Coulomb correlation in presence of antiferromagnetic 
(AFM) coupling is responsible for the insulating ground state. In the ground state the d xz and d yz 
orbitals are ordered and intra-chain vanadium ions are antiferromagnetically coupled. The calcula- 
tion gives small spin-orbit coupling (SOC), which provides a tilt of ~11.3° to the magnetic moment 
from the z-axis. In the presence of weak SOC and strong exchange coupling, the experimentally 
observed small magnetic moment and low AFM transition temperature appear to arise from spin 
fluctuation due to activeness of geometrical frustration. 

PACS numbers: 75.25.Dk, 71.20.-b, 71.27.+a 



I. INTRODUCTION 

The orbital degree of freedom is an important entity in 
the condensed matter physics which plays a crucial role 
in stabilizing many exotic phases observed in the strongly 
correlated electron systems. 1 When the degenerate d or- 
bitals of the transition metals are partially filled then 
occupation of a particular orbital at a particular site is 
expected to be dictated by the occupation of another or- 
bital at another site, which can lead to various kind of 
orbital ordering (00) similar to the spin ordering. For 
example, in the case of LaMnOs degenerate e g orbitals 
are occupied by only one electron and predicted to show 
antiferro-orbital ordering (AFOO) in the xy plane in 1955 
by Goodenough. 2 The Coulomb correlation is found to 
play an important role in the orbital physics of transi- 
tion metal oxides. However, it is still not clear whether 
it just enhances the effect of lattice distortions or really 
drives the 00 via superexchange^ In spite of the am- 
biguity about the exact mechanism behind the 00, it is 
almost clear that such 00 is often accompanied by reduc- 
tion in the crystal symmetry. Thus in the geometrically 
frustrated system 00 is expected to relax the frustration 
leading to the formation of novel magnetic phases earlier 
forbidden by the frustration. 

Spinel vanadates with general formula AV2O4 (A-Cd, 
Mg, and Zn) is an important geometrically frustrated 
system which has attracted a great deal of attention for 
a decade because of 00 induced structural transition 
and formation of fascinating magnetic phases A"— All the 
studied compounds show cubic to tetragonal transition at 
low temperature and paramagnetic (PM) to AFM tran- 
sition at slightly lower temperature^"—^ The consensus 
is emerging among the researchers regarding the 00 in- 
duced structural transition, however the exact pattern 
of 00 remains a matter of controversy ! WW 5 More- 
over, the exact role of SOC along with its strength in 
stabilizing the magnetic and orbital ordering is yet to be 
decided in this series of compounds A 1Q i 14 In the present 
work we would like to address these issues for MgV2 04 



(MVO). The recent experimental work on this compound 
has shown quite different results in comparison to the 
well studied ZnV2 04 (ZVO). 15 In the tetragonal phase 
the space group of MVO is I4m2 whereas that of ZVO 
is I4i/amd. The magnetic moment (MM) of V ion in 
MVO is ~0.47 fiB which is ~0.15 /i# less than that in 
ZVO. Such a reduced value of MM in ZVO is due to large 
but negative contribution from the orbital part of MM. 14 
However, experimental data of MVO do not suggest such 
a large contribution from the orbital part and indicating 
to a deeper reason for observed small MM. Moreover, the 
MM is also seen to make an angle of ~8° with the z-axis 
indicating to a weak SOC in MVO compound. 

Here we explore the role played by spin and orbital de- 
grees of freedom in deciding the electronic and magnetic 
properties of MVO by using ab initio electronic struc- 
ture calculations. The AFM interaction in the presence 
of strong Coulomb correlation is found to be crucial in 
driving the system to insulating ground state. The d xz 
and d yz orbitals get ordered in the tetragonal phase and 
00 becomes more robust in the presence of AFM in- 
teraction. The spin and orbital part of MM is found to 
~1.4 and -0.2 /ig, respectively. In the light of this re- 
sult and in the presence of large exchange coupling of 
~58 meV one can suggest that the geometrical frustra- 
tion may be responsible for the experimentally observed 
low AFM transition temperature (~42 K) and small MM 
(-0.47 fji B ). 



II. COMPUTATIONAL DETAILS 

The non-magnetic (NM), ferromagnetic (FM), and 
AFM solutions of MVO are obtained by using state-of- 
the-art full-potential linearized augmented plane wave 
(FP-LAPW) method. 16 The lattice parameters and 
atomic positions used in the calculations are taken from 
the literature^ The muffin-tin sphere radii automati- 
cally set in the calculations are 1.5664, 1.7045, and 1.4943 
Bohr for Mg, V, and atoms, respectively. For the ex- 
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change correlation functional, we have adopted recently 
developed generalized gradient approximation (GGA) 
form of Perdew et alJ^ The effect of on-site Coulomb in- 
teraction is also considered within GGA+/7 formulation 
of the density functional theory. 18 In GGA+/7 method 
the U and J are used as parameters. We varied U from 
3-5 eV and fixed J=0.5 eV. We found similar results for 
all values of U. The results correspond to £7=4 eV are 
only discussed in the manuscript. In order to see the role 
of orbital degrees of freedom on the electronic and mag- 
netic properties of the compound SOC is also considered 
in the calculations. The self-consistency was achieved 
by demanding the convergence of the total energy to be 
smaller than 10 -4 Hartree/cell. 



III. RESULTS AND DISCUSSIONS 

The atomic arrangements in the unit cell are shown in 
Fig. 1. It is evident from Fig. 1(a) that each V atom 
is surrounded by six O atoms forming an octahedron. 
The octahedra are edge shared to each other. The small 
trigonal distortion splits six V-0 bonds of an octahedron 
in two groups containing three bonds each with bond 
length of 2.016 and 2.033 A. The four nearby V atoms 
form a regular tetrahedron with an edge of length 2.98 A. 
Each of the tetrahedron is surrounded by four neighbor- 
ing tetrahedra via corner sharing and forming the chains 
of V atoms, see Fig. 1(b). In the tetragonal phase the 
tetrahedra become distorted with the V-V bond length 
of 2.971 and 2.980 A. Such a small distortion would re- 
duce the geometrical frustration and is expected to give 
rise to novel electronic and magnetic phases. 

In order to know the exact ground state of the com- 
pound we obtained various magnetic solutions using 
tetragonal structure. The total density of states (TDOS) 
correspond to these solutions are shown in Fig. 2. It is 
clear from Figs. 2(a)-(d) that the GGA solutions provide 
metallic state as opposed experimentally observed insu- 
lating behavior. This result is not surprising as the GGA 
underestimates the Coulomb correlation among the 3d 
electrons which is often found to be responsible for in- 
sulating ground state of the transition metal oxides J^r— 
The NM and FM solutions obtained from GGA+Z7 cal- 
culations also provide metallic state as evident from Figs. 
2(e) and (f). This indicates that there may be a deeper 
reason for the insulating ground state of the system. At 
this juncture, it is important to note that the FM solu- 
tion of ZVO within LSDA+Z7 is found to derive insulat- 
ing ground state of the compound*^ In order to know 
the exact cause for insulating ground state of MVO we 
performed AFM calculations. The AFM solution creates 
a soft gap and in presence of SOC it provides a hard gap 
of —0.16 eV, see Figs. 2(g) and (h). It is important to 
note that the increased value of U enhances the band 
gap in both the cases. Thus present work clearly estab- 
lishes that the AFM coupling of V moments in presence 
of strong on-site Coulomb interaction is responsible for 



insulating ground state and SOC provides robustness to 
the insulating property of the compound. 

Now we discuss the effect of different interaction pa- 
rameters on the electronic structure of the compound. 
First we start with GGA results. In the absence of mag- 
netic interaction there are large density of states (DOS) 
at the Fermi level (ep) with dominating contribution 
from the V 3d states. This may be considered as a signa- 
ture of magnetic ground state under Stoner theory. The 
magnetic interaction reduces the DOS at ep by —50% due 
to exchange splitting of the bands contributing at the e p 
and providing almost half- metallic state, see Fig. 2(b). 
Moreover, the energy of FM solution is found to be —0.65 
eV/fu (fu=formula unit) less than that of NM solution 
indicating the magnetic ground state. The energy differ- 
ence between the band edge of the up and down spins 
may be considered as a measure of exchange interaction 
which is found to be -0.4 eV. The MM of V is found to 
be —1.25 hb- Interestingly, FM interaction induces finite 
MM (-0.14 fj, B ) at the Mg ions occupying 2c (0,1/2,1/4) 
Wyckoff sites. The total MM/fu comes out to be -3.6 
which corresponds to 5~1 state of the V 3+ ion. It is 
evident from Fig. 2(c) that the AFM interaction among 
the V moments reduces the DOS at the ep drastically 
(—4 times less than that of FM). The ep lies at the min- 
ima of DOS which is a reminiscence of the pseudo-gap. 
The AFM interaction reduces the band width (BW) of 
the system and the BW of the deeper bands decreases by 
-0.3 eV. Further, it decreases the MM of V by -0.1 \i B 
and does not create any MM at Mg sites. The energy /fu 
of the AFM solution is -0.16 eV less than that of FM so- 
lution. This is a clear evidence of the AFM ground state 
whose spin ordering will be discussed in the later part of 
the manuscript. The inclusion of SOC at this stage does 
not have any significant effect on the electronic structure 
of the compound as evident from Fig. 2(d). By com- 
paring the energy of the AFM and AFM-j-SOC solutions 
one can get the rough estimate of SOC strength of V 3d 
electron as the contribution from Mg and O atoms is ex- 
pected to be negligibly small. Our GGA calculation gives 
the strength of SOC of —5 meV for the V 3d electrons. 

Figs. 2(e)-(h) depict the effect of on-site Coulomb 
interaction among the V 3d electrons on the electronic 
properties of the compound in the presence of various 
interaction parameters. On-site Coulomb correlation re- 
duces the BW as it localizes the electrons. In the case of 
NM and FM solutions there is a drastic decrease in the 
V 3d DOS at the ep due to transfer of spectral weight 
(earlier contributing at ep) away from it. Moreover, FM 
interaction gives rise to a perfectly half-metallic state 
with a band gap of —2.5 eV in the down-spin channel. 
In the presence of AFM interaction among the V spins 
the system becomes insulating due to formation of upper 
and lower Hubbard bands. Further, inclusion of SOC in- 
creases the separation between upper and lower Hubbard 
bands. The exchange interaction estimated from the FM 
solution is —0.5 eV which is 0.1 eV more than that ob- 
tained from simple GGA calculation. This enhancement 
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is attributed to increased Hund's coupling strength due 
to increased spatial localization of V 3d electrons, which 
also enhances the MM at V sites by ~0.2 \ib- The to- 
tal MM/fu comes out to be 4 which corresponds to 
5=1 state of the V 3+ ion. Such a large value of mag- 
netic moment at V site is in sharp contradiction with the 
experimentally observed small magnetic moment. 15 This 
clearly indicates that some other parameters are playing 
important role in deciding the magnetic properties of the 
compound. On comparing the energy of various solutions 
we find that the AFM state is the true ground state of 
the system as energy of NM>FM> AFM. 

In order to study the role of orbital degrees of freedom 
we have performed FM GGA+Z7 calculations in both cu- 
bic and tetragonal phases. In the cubic phase the occu- 
pancies of d xz and d yz orbitals are same at every V sites, 
whereas for the tetragonal phase the occupancies of these 
orbitals are found to be different at different site which 
is a direct evidence of 00 taking place in the tetragonal 
structure. Moreover, the 00 pattern does not depend 
on the nature of magnetic interaction as evident from 
Table 1 where we have listed the occupancy of d x 2_ y 2^ 
d XZl and d yz orbitals of four V atoms forming the tetra- 
hedron and obtained from FM and AFM solutions. It 
is evident from the table that each site is occupied by 
d x 2_ y 2 orbital. The VI and V2 sites are mainly occupied 
by d xz orbital and that of V3 and V4 sites by d yz . At 
this juncture it is important to note that 00 is observed 
in the PM phase of the spinel vanadates and in the PM 
phase there is a local MM at the V site and hence mag- 
netic solutions would provide the better representation 
of the PM state in comparison to NM solution. Table 
1 also indicates the AFM coupling between VI and V2 
(V3 and V4) and FM coupling between VI and V4 (V2 
and V3). Interestingly, AFM interaction appears to pro- 
vide more stability to the 00 as the occupancy of d yz 
{dxz) orbital at VI and V2 (V3 and V4) sites is found 
to decrease by ~0.08. Moreover the energy/fu of AFM 
solution is also ~0.24 eV less than that of FM solution 
suggesting the AFM ground state. This energy difference 
between AFM and FM solutions is ~80 meV less than 
that obtained from GGA solution. This highlights the 
importance of Coulomb correlation in establishing the 
AFM ground state. The MM of vanadium ions corre- 
sponds to FM and AFM solutions are found to be ~1.43 
and 1.35 respectively. 

As mentioned in the introduction that the 00 is nor- 
mally considered as a cause for the structural transition 
in the PM phase of spinel vanadates. However, based 
on present work it is difficult to say whether 00 is the 
cause of structural transition or it is just an effect of it. 
In order to understand the cause of 00 seen in different 
transition metal oxides mainly two mechanism exist in 
the literature which are purely electronic and structural 
in origin* 1 -^ 4 - To separate out these two contributions to 
the 00, Pavarini et al. have carried out beautiful work 
on two canonical 00 systems viz. KCUF3 and LaMnOs, 
where they have used LDA+DMFT method^ 4 It is im- 



portant to note that GGA+?7 method used in the present 
work is a static mean-field theory whereas LDA+DMFT 
used in the work of Pavarini et al. is a dynamical mean- 
field theory and hence better in approximation. Thus in 
order to know the exact cause of 00 in the MVO com- 
pound work in line with Pavarini et al. is desirable. 

Now we study the effect of SOC on the magnetic state 
of the compound. The GGA+Z7 +SOC solutions also give 
AFM ground state as the energy/fu of AFM solution is 
found to be ~0.18 eV less than that of FM solution. As 
mention above the AFM interaction provides more stabil- 
ity to the 00 which further enhances the orbital moment 
as evident from Table 2 where we have shown the spin 
(5), orbital (L) and total (J) moment of the V ion cor- 
responds to FM and AFM solutions. The orbital part of 
MM in the AFM state comes out to be ~-0.2 /ib which 
is ~7 times less than the spin part of MM suggesting the 
weak SOC in MVO with respect to ZVO 14 where large 
orbital moment of -0.75 \ib has been reported. The di- 
rection of total MM is found to be ~11.3° away from 
the z-axis. The small value of orbital moment and the 
direction of total magnetic moment are in consonance 
with the experimental findings where neutron scattering 
studies have revealed the small orbital moment and MM 
is tilted at ~8° from the z-axisji 5 - However, the above 
calculated values of spin and orbital moments cannot ac- 
count for the experimentally observed magnetic moment 
of 0.47 jiB- At this point it is important to note that the 
experimentally estimated magnetic moment of 0.63 \±b 
for ZVO is well accounted by taking into account the cal- 
culated large but negative orbital moment of -0.75 /i£. 14 
This indicates that the orbital sector of the MVO is not 
as influential as found in the ZVO in deciding the mag- 
netic state of the V ions. Thus there may be a deeper 
reason for the experimentally observed small magnetic 
moment in the MVO which will be discussed in the next 
paragraph. The final spin and orbital ordering patterns 
obtained from the calculations are shown in Fig. 3. The 
spins are forming AFM chains in the x and y directions 
and nearest neighbor AFM chains are connected by lines 
with FM ordering. The AFM chains are accompanied by 
ferro-orbital ordering (F00) where d xz and d yz orbitals 
are occupied along the x and y directions, respectively. 
The AFOO supports the formation of FM chains where 
neighboring sites are alternatively occupied by d xz and 
d yz orbitals. These spin and orbital ordering patterns are 
in accordance with the Goodenough-Kanamori schemes. 
Here it is important to note that historically Goodenough 
has given the semi-covalent scheme for explaining the 
experimentally observed complex magnetic structures in 
Lai-^Ca^MnOa and predicted different 00 corresponds 
to different spin arrangements. 2 According to this scheme 
length of FM bond should be greater than that of AFM 
bond. However, we have observed opposite behavior as 
AFM bond is 0.03 A larger than the FM bonds. 

As mentioned above that the calculated small orbital 
moment of -0.2 \±b for MVO cannot account for the 
experimentally observed value of total MM ^0.47 
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whereas calculated large orbital moment of -0.75 hb for 
ZVO provides a good description of its experimentally 
observed total MM of 0.63 [lb> These results appear 
to suggest that the geometrical frustration is still active 
in the tetragonal phase of MVO which can give rise to 
spin fluctuation at low temperature. Such spin fluctua- 
tion is expected to reduce the MM drastically. The level 
of frustration in magnetic systems is defined by frustra- 
tion index / = \9cw\/T*, where Ocw is the Curie- Weiss 
temperature and T* is the critical temperature at which 
the system ultimately develops long-range spin order. 24 
Higher the value of / more will be the level of frustra- 
tion. Thus the above conjecture about the activeness 
of frustration can further be tested by estimating the 
Heisenberg exchange interaction strength (Jh) between 
V moments on which Ocw depends. The rough estimate 
of it can be found by mapping the energies of FM and 
AFM solutions to the Heisenberg Hamiltonian. 25 Our cal- 
culation gives Jh~58 meV. Using this value of Jh we 
have estimated the AFM transition temperature based on 
mean- field theory and it comes out to be ~925 K, which 
is closer to the experimental value of 6 C w^ Using exper- 
imental T*^42 K we have estimated the value of /«22. 



Such a large value of / further strengthens the conjecture 
about the activeness of frustration in MVO compound. 
Here it is important to note that the GGA+Z7 formula- 
tion of density functional theory is a mean-field theory, 
which is not capable of addressing issue related with spin 
fluctuations directly. Thus, to address this one needs to 
go beyond the mean-field theory. 

IV. CONCLUSIONS 

In conclusion, we have investigated the electronic 
and magnetic properties of a geometrically frustrated 
MgV204 by using ab initio electronic structure calcu- 
lations. This compound is a Mott-insulator and its insu- 
lating ground state is arising due to combined effect of 
strong Coulomb correlation and AFM interaction. The 
d xz and d yz orbitals are found to be ordered in the tetrag- 
onal phase. The spins are forming AFM ordered chains 
along the x and y directions and making an angle ~11.3° 
with the z-axis. The SOC is weak and geometrical frus- 
tration appears to be active in deciding the magnetic 
state of the system. 
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FIG. 1. (Color online) Atomic arrangement of the unit 
cell. The formation of V06 octahedra and V4 tetrahedra 
are shown in (a) and (b), respectively. 

FIG. 2. (Color online) Evolution of total density 
of states (TDOS) with various interaction parameters. 
Please see the text for the details. 

FIG. 3. (Color online) (a) Intra-chain antiferromag- 
netic ordering along the x and y directions, (b) spin and 
orbital arrangements at the tetrahedron level. 

Table 1: Occupancies of d x 2_ y 2, d xz , and d yz orbitals 
and magnetic moments of four V atoms (viz. VI, V2, 
V3, and V4 forming the tetrahedron, see Fig. 3(b)) cor- 
responding to FM and AFM (in brackets) solutions ob- 
tained from GGA+U (E/=4 eV) calculations. 
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VI 


V2 


V3 


V4 


dx 2 —y 2 


0.62(0.65) 


0.62(0.65) 


62(0.65) 


62(0.65) 




0.64(0.66) 


0.64(0.66) 


0.18(0.1) 


0.18(0.1) 


dy Z 


0.18(0.1) 


0.18(0.1) 


0.64(0.66) 


0.64(0.66) 


MM (fi B ) 


1.43(1.35) 


1.43(-1.35) 


1.43(-1.35) 


1.43(1.35) 



Table 2: The expectation value y and z components 
of spin (£), orbital (L) and total (J) moment of V ion ob- 
tained from FM and AFM (in brackets) GGA+/7+SOC 
(Z7=4 eV) solutions. 
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-0(-0) 


-0(-0) 


0.71(0.67) 


L 


-0.02(-0.1) 


-0(-0) 


-0.07(-0.17) 


J 


-0.02(-0.1) 


-0(-0) 


0.64(0.5) 
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